Molecular theory of elastic constants of liquid crystals. III. 
Application to smectic phases with tilted orientational order 

Yashwant Singh and Jokhan Ram 

Department of Physics, Banaras Hindu University, 
Varanasi 221 005, India 
(February 1, 2008) 



Abstract 



Using the density functional formalism we derive expression for the distor- 
tion free energy for systems with continuous broken symmetry and use it to 
derive expression for the elastic constants of smectic phases in which director 
is tilted with respect to the smectic layer normal. As in the previous papers of 
the series ( Phys. Rev. A 45, 974 (1992), E 49, 501, (1994) ) the expressions 
for the elastic constants are written in terms of order and structural param- 
eters. The structural parameters involve the generalised spherical harmonic 
coefficients of the direct pair correlation function of an effective isotropic liq- 
uid. The density of this effective isotropic liquid depends on the nature and 
amount of ordering present in the system and is evaluated self- consistently. 
We estimate the value of elastic constants using reasonable guess for the order 

and structural- parameters. 
61.30 Cz, 62.20 Di, 61.30Jf 
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I. INTRODUCTION 



In a previous paper |IJ of this series ( hereafter referred to as I ) we developed a theory 
based on the density functional formalism for the deformation free energy of any system 
with continuous broken symmetry and applied the theory to derive expressions for the 
elastic constants of uniaxial nematic and smectic A (Sm A) phases. In the second paper 
H (hereafter is referred to as II) the theory was applied to derive expression for all the 
twelve bulk elastic constants of biaxial nematic phase of orthorhombic symmetry. These 
expressions of elastic constants are written in terms of order parameters which characterize 
the nature and amount of ordering in the phase and the structural parameters which involve 
the generalized spherical harmonic coefficients of the direct pair correlation function of an 
effective isotropic liquid the density of which is determined using a criterion of the weighted 
density-functional formalism. The purpose of this paper is to expand upon the theory given 
in I and apply it to those smectic phases in which the director is tilted with respect to the 
smectic layer normal. This class of smectics include Sm C, Sm I, Sm F and Sm K phase and 
their chiral versions. 

In general, smectic liquid crystals have a stratified structure with the long axes of the 
rod-like or lath-like molecules parallel to each other in the layers. This situation corresponds 
to partial breakdown of translational invariance. Since a variety of molecular arrangements 
are possible within each layer a number of smectic phases are possible. Thus a smectic is 
defined by its periodicity in one direction of space and by its point group symmetry which 
is continuous or discrete subgroup of -Dooh. In this article we are confined to those smectic 
phases which point group symmetries are C^h and Ci- Most important phases of this class 
of smectics as far as elasticity is concerned, are Sm C and SmC* phases. 

In Sm C the molecules are arranged as shown in Fig.l. Each smectic layer is a two 
dimensional liquid. The director n makes a finite angle with the normal to the layer. The 
tilted orientational order is described by a two dimensional unit vector which is along the 
projection of n on the layer and is denoted as c ( see Fig.l). The direction of c is defined 
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relative to a chosen axis in the plane and simultaneously with the sign of projection of n 
on the normal to the layer. The correlation of the tilt direction of different layers implies 
that c must be uniform over macroscopic distance. Though Sm C is n <-> — n symmetric, 
such symmetry operation does not hold for the c -director. The states described by c and 
— c are not equivalent and hence c is a two dimensional polar vector. The tilted structure 
introduces biaxiallity within the layers and Sm C is optically biaxial with point symmetry 

The other non-chiral smectic phases of this class namely Sm I, Sm F and Sm K all 
possess more order than the Sm C phase. In these each smectic layer is an independent 
two dimensional bond orientational ordered system. The coupling between neighbouring 
hexatic layers drives the quasi long range hexatic bond-orientational ordering of a layer into 
a truly long range ordered hexatic state. Therefore these phases have long ranged bond- 
orientational order but short ranged positional order in the layer in a three dimensional 
stacked hexatic systems. The director n, as in Sm C, makes a finite angle with the layer 
normal. These phases, in spite of having quite different local order, have the same symmetry 
as a Sm C phase. The phase transition between Sm C and one of the tilted hexatic phases 
is found to be first order. This transition is analogous to the 3-dimensional liquid to vapor 
transition which is first order upto a critical point, and beyond the critical point there is no 
real phase transition. One can go continuously from one to the other phase without phase 
transition at all. The expression for the elastic free energy density for these phases are, 
therefore, similar to that of Sm C; the difference can be only quantitative. In view of this 
we discuss the elasticity of Sm C phase only. 

The chiral smectic C (SmC*) has modulated structure at a scale dimension of the order 
of 1 ^m and larger. The modulated (helical) structure occurs as a result of a precession 
of the molecular tilt about an axis perpendicular to the layer planes. The tilt direction is 
rotated through an azimuthal angle on moving from one layer to the next. As the rotation 
is in a constant direction, a helix is formed which is either left or right handed. In SmC* the 
mirror plane is lost but C2 axis is present locally in each smectic layer (it is perpendicular to 

3 



the plane containing the molecular axes). The C 2 axis in the case of SmC* phase is a polar 
axis which admits the existence of a spontaneous polarization along it. The SmC* phase 
therefore exhibits ferroelectricity. 

Many compounds exhibit phases such as antiferroelectric (SmC^), ferroelectric (Sm C*) 
etc in addition to SmC*, as the temperature is varied. The director n in the SmC^ phase, 
like in SmC* phase is uniformly tilted with respect to the layer normal. However, unlike in 
SmC* phase, the difference in azimuthal angle of n in successive layer is ir + a. The small 
angle a gives rise to helical structure of n along the layer normal and arises from the chirality 
of the constituent molecules. In SmC* the layers are stacked in such a way that there is a 
net overall (in contrast to SmC*^ phase) spontaneous polarization. This means the number 
of layers of opposite polarization is not equal. The elasticity associated with the structure 
of all chiral smectic phases is same as that of SmC* phase. 

The paper is organised as follows : In Sec. II we describe the elastic continuum theory 
for Sm C and SmC* phases and derive expression for the elastic free energy density. This 
expression is then written in terms of deformation variables so that it becomes appropriate 
for comparison with an expression found from microscopic theory. The density functional 
approach has been used in sec. Ill to describe the deformed state of smectics with tilted 
orientational order. The expression for the deformation free energy density derived here is 
more general than that given in I in the sense that all the cases discussed in I correspond to 
tilt angle equal to zero. The expressions for the elastic constants are derived by comparing 
the elastic free energy found from the density functional theory with that given in Sec II. 
These expressions of the elastic constants involve the order parameters which measure the 
nature and amount of ordering in the system and direct pair correlation function (DPCF) of 
an effective isotropic liquid. In Sec. IV we discuss the relative magnitude of these constants 
using reasonable guess for the order parameters and the spherical harmonic coefficients of 
the DPCF. 
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II. CONTINUUM THEORY FOR ELASTIC FREE ENERGY DENSITY 



An undeformed smectic has parallel and equidistant layers and position independent 
director. When some small amplitude long wavelength distortions are imposed on this ideal 
state, layers may get displaced and curved. Let u(x,y,z) represent displacements of the 
layer normal to their planes. A layer at z before displacement is now at 

z(x,y) = z + u(x,y, z) (2.1) 

We choose a space-fixed (SF) frame S° such that its z-axis is along the normal to the 
unperturbed smectic layers (see Fig.l). The system is described by (i) the angle Q z of 
director c with respect to Ox-axis and (ii) the vertical displacement u of the layers whose 
derivatives 

du du 
Q x = — and Q y = -— (2.2) 

represent small angles of rotation about x and y-axes, respectively. 

Since the free energy of a system with broken continuous symmetry is invariant with 
respect to spatially uniform displacements and rotations that take the system from one 
point in the ground state manifold to another, the elastic free energy must be function of 
derivatives of Q or second order derivatives of u which correspond to the curvature of the 
planes. The contribution to free energy density due to deformation is called elastic free 
energy density / e (r). To find f e (r) one expands f{r) that belongs to the deformed state 
around the ideal state having the free energy density f u . The expansion is expressed in 
terms of the spatial derivative of the order parameter fields. The elastic continuum theory 
deals only with small spatial derivatives. Consequently, only the lowest order terms in the 
expansions are taken into account. 

Because of the symmetry of Sm C the free energy density must be invariant of the 
simultaneous transformations x — ► —x, y — > — y and z — > —z (due to centre of symmetry) 
and by the transformation y — > —y as the vertical xOz plane is the plane of symmetry. Thus 
for the elastic free energy density one has 



fe(r) = fdc{r) + f dk (r) + f dck {r) + f dl (r) (2.3) 

where f dc represents contribution arising due to distortion of director c while f dk is associated 
with curvature of layers. j dck represents the contribution arising due to coupling of different 
deformation modes whereas f d i represents the contribution due to variation in layer thickness 
as in Sm A. Thus |[J 

hi = \Bl 2 (2.4) 
where 7 = describes the variation in the thickness of layers. 

To derive expressions for f dc , f dk , f dck we consider director triad formed by three orthonor- 
mal unit vectors c, k, p and let at origin the director c is along the x-axis of SF frame, k 
along z - axis, i. e. along normal to the layers and p along y-axis, i. e. 

c = (1,0,0), Po = (0,1,0), k = (0,0,1) (2.5) 

The orientation of the director triad at a neighbouring point is rotated with respect to 
the space-fixed coordinate frame due to deformation and is given as 

c=(l,c y ,c z ), p = (p x ,l,p z ), k=(k x ,k y ,l) (2.6) 

Because of orthonormality condition of director triad 

Px = ~ c yi Pz = ~kyi c z = ~kx (2-7) 

For small deviations we have 

dy v v dx v y ' 

Combining the three orthonormal vectors c, k, p we can form nine invariants. The 
elastic free energy density can be expressed in terms of these invariants || 

dVt x 

Du = Cappdakp = c.Vk.p = — — 
D 2 i = PaPpdakp = p.Vk.p = 



An 


= k a Ppd a kp = 


k.Vk.p = - 


Dn 


1 O 

= c a kpd a cp = 


-c.Vk.c = 


B> 22 


^ U £1 n 


— p.VK.C = 


D 32 


= kakpdaCp = 


-k.Vk.c = 


D 13 


= C a Cpd a P{3 = 


— c.Vc.p = 


D 23 


= p a cpd a p p = 


— p.Vc.p = 


D 33 


= k a cpd a p p = 


-k.Vc.p = 



dVt x d 2 u 



dz dxdz 

dVty 



dx 

dily _ 

dy 

dQ„ d 2 u 



dz dydz 

on, 



dx 
dy 

on, 



(2.9) 



dz 

Note that the invariants D 3 i, and D 32 correspond to variation of layer compression (or 
dilation) along y and x-axes, respectively : 

D 31 = 9Qx = — = -?1 = -p V7 
dz dydz dy 

D 32 = _«k = -£L = _»T = -c.V 7 (2.10) 

dz dxdz dx 

These terms are omitted as their contributions are small compared to linear term f dt (see 
Eqn. (2.4)). We are therefore left with six independent invariants. 

The elastic free energy density due to orientational deformation of Sm C can be written in 
terms of these six independent invariants. Each term should, however, satisfy the symmetry 
elements of Sm C. Note that (k, c) plane is a symmetry plane, p is C 2 axis and there is also 
inversion symmetry. We can therefore write elastic free energy density as a sum of square of 
each invariant and their products taken in such a way that k and c enter the expansion in 
an even number of times (i.e. k and c are considered without discrimination) and p which 
also appear even number of times but on its own. Thus 

f dc = \b x D\ z + \b 2 D\ z + \b 3 D\ z + B 13 D 13 D 33 (2.11a) 

f dk = \a(D 11 + D 22 f + l -A l2 D\ 2 + l -A 2l D 2 2l (2.11b) 

Uk = C^uDn + C 2 D 21 D 23 + C 3 D 12 D 23 + C 4 D 12 D 21 (2.11c) 

Since 



Dl2 D 23 = (*?-) = » („.«?.) _ » ( n ,*V) + «k«k, ( 2.i2) 

\ ox J \ oy J oy \ ox J ox \ oy J ox oy 

it is the sum of two gradient terms (surface terms) and a term analogous to D n D 13 . Similar 
decomposition can also be made for D X2 D 21 . We therefore have only two independent terms 
in the expression of fdck due to coupling between the deformation of director c and the 
curvature of smectic layers, as given in ref. [4]. 

Since fa c involves director c and can be written in terms of (V.c) 2 , (c.V x c) 2 and 
(c x V x c) 2 as 

(v ' d)2 =(f) 2 = ^ 

and 

(k.V x £)(c.V x £) = = D 13 D 33 (2.13) 

Therefore the constants B 1: B 2 , B 3 , B 13 can be identified as elastic constants associated 
with bend, splay, twist and coupled twist-bend modes of deformations of director c and 
have the dimension of energy per unit length (dynes). The last term in Eq.(2.11a) is allowed 
because of the Sm C is invarient under the simulataneous transformation n — > — n and 
c — > — c. These constants are expected to be comparable in magnitude to the Frank elastic 
constants of a homologous nematic phase. 

The contribution f dk is analogous to Sm A and is due to layer deformation (curvature). 
The term A, A\ 2 and A 2 \ represent splay contribution of a system of a monoclinic symmetry 
as Sm C has C 2 h symmetry in contrast to of Sm A. The value of these constants 

are therefore expected to be of the same order of magnitude as in Sm A phase. B (see 
Eq.(2.4)) which has dimension of energy per unit volume is associated with possible change 
of interlayer distance. 

Since Sm C* has no inversion symmetry, its elastic free energy density may have terms 
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linear in invariants of Eq (2.9). However, due to symmetry requirements only those invariants 
which involve k and c in pairs without discrimination and p alone contribute to linear terms. 
Thus 

jf ] = -£>!(£. Vc.p) + D 2 (p.Vk.c) - £> 3 (k.Vc.p) (2.14) 
where superscript (1) stands for linear terms. In terms of f2 vector field 

Out of the three terms in Eq.(2.15) the term D 3 is most important as it corresponds to 
simple twist of director c. The terms D\ corresponds to bend of director c while D 2 tends 
to transform a flat layer into a twisted ribbon. The terms D 1 and D 2 are expected to be 
negligible in comparison to D 3 term. 

If g = — (-§^) where B 3 is defined in Eq.(2.11a) the contribution fd c (r) for Sm C* is 
modified to 

« r) = 2 B '(afj + 2 B 'i^) + 2*(^ + ?0 ) +Bl3 l^)(lfc) (216) 

while other terms of Eq.(2.11) remain same. The ground state of Sm C* has a helical 
structure. The B% term can therefore be considered to be a contribution due to deformation 
from this uniform helical structure. 

The expressions for the elastic constants are obtained by comparing the expression of 
elastic free energy density obtained from the continuum theory and from a microscopic 
theory involving molecular parameters. In order to make this comparison possible we have 
to rewrite Eq.(2.16) in terms of distortion vectors defined below. 

Since in the ground state the directors are position independent we take S* to be the 
unperturbed director frame at a point R. Let S* be described by the unit vectors such 
that (see Eq. (2.5)) 

N« = c , N 2) = p , N( 3) =k , (2.17) 
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Let S be the perturbed director frame described by the unit vector such that (see 
Eq. (2.6)) 

NW = c = c + (t( 1 ), N( 2 ) = p = p + ^ 2 ), N( 3 ) = k = k + a( 3 ), (2.18) 

where may be regarded as the distortion vectors, and let 

N^N«.<» = ^, + ( 7« (2.19) 

be the components of the perturbed directors with respect to S° frame. 

The transformation from S* to S is effected by 3 x 3 orthogonal matrix T with element 

T jk = Nj?> = 5 jk + o? 

such that 

= 7} fe Nf = N?' } + o^Nf 
The orthogonality condition TijT ik = 5j k leads to 

of = -<r ( £ ) -vf^ for j ^ k 



for i = * ( 2 - 2 °) 

where prime on summation indicates no summation over j. Thus for small deformations 



of ~ O and of ~ O 2 

where O means first order small quantity. 

From above relations it follows that of, o^ 3 " 1 and of can be regarded as basic indepen- 
dent distortion variables; all as a function of R. Every other component can be constructed 
with the help of Eq.(2.20). 

From above discussions we find 

k x = — ^1 

k y = -Q x = -of (2.21) 
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Using these relations we rewrite Eq. (2.11) as 
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Fdc = I dK 



a)' 



a)' 



+ B 



13 



da?' 
dx 



day 

dz 



Fdk — I dH 



dck 



1 A da± 
2 I dx 



(2)\ 2 1 /a_(3) x 2 



l A da? 
2 12 ~dx~ 



1 a / ^3 



(2)' 



/ 



dR 



Ci 



9a; 



+ C 2 



(2.22a) 

(2.22b) 
(2.22c) 



9y I \ dy 

It is often convenient to assume long wavelength deformations in space fixed frame having 
the form 



<7?> = Bf sin(q.R) 



(2.23) 



where B*jp are small amplitudes and q is the wave vector. Substituting this into Eq. (2.22) 
leads to 



Fdc 

~v 

Fdk 

~v 

Fdck 
V 



\ [bM ] )V x + B 2 (Bi 1] )V y + B,{B?fql + B^fq^ 

^[^ 2) ) 2 ^ + ^2(4 3) ) 2 ^+A 2 l(4 2) )^ 

\c^B?)ql + \c 2 {B^B^l 



(2.24a) 
(2.24b) 
(2.24c) 



In writing above equations we have used the relations 

J dRsin 2 (q.R) = J dRcos 2 (q.R) = 



V 
~2 



and 



J dRsin(q.R)cos(q.R) = 



III. DENSITY-FUNCTIONAL APPROACH 
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A. Expression for the distortion free energy 



A liquid crystalline phase is characterized by the spatial and orientational configuration 
of molecules. At the phase transition point these configurations undergo a modification, i.e., 
abrupt change may take place in the symmetries of the system. The molecular configurations 
of most ordered phases are adequately described by the single particle density distribution 
p(x). The vector x is taken here to indicate both the location r of the centre of a molecule 
and its relative orientation f2 described by Euler angles 0, 9, ip. For an isotropic uniform 
system p(x) is independent of positions and orientations. 

The single particle density distribution p(x) provides us with a convenient variational 
quantity to specify an arbitrary state of a system. One may consider a variational thermo- 
dynamic potential W(T, P, [p(x)]) as a function of p(x). The equilibrium state of the system 
at a given T and P is described by the density p(T, P, x) corresponding to the minimum of 
W with respect to p(x). This forms the basis of the density- functional theory 0. 

The basic thermodynamic potential used to determine the isothermal elastic properties 
of a system consisting of N particles in volume V at temperature T, is the Helmholtz 
free energy A[p\. Elasticity is associated with the behaviour of A[p] with respect to small 
deformation of the system away from its equilibrium (ground) state [|TJ. 

In the density-functional formalism the free energy of a system is expressed in terms of 
the direct correlation function of the medium ||. 



and /3AA[p] is the excess reduced free energy arising due to intermolecular interactions. 
Here j3 = (fc^T) -1 where ks is the Boltzmann constant and T is the temperature. In the 
weighted density-functional approach 



PA\p]=l3A id \p]+/3AA[p] 



(3.1) 



where the ideal gas part 




(3.2) 



12 



/3AA[p] = -- J tfoti J dx 2 p(xi)c(x 1 ,x 2 ,p)p(x 2 ) (3.3) 

Here the function c represents the direct pair correlation function (DPCF) of an effective 
isotropic fluid. The effective density p is found using the relation || 

p[p] = J rfxi J dx 2 p(xi) p(x 2 ) ^(xi,x 2 ;p) (3.4) 

where po is the averaged density of the ordered phase and u is a weight factor. p[p] is viewed 
here as a functional of p(x). To ensure that the approximation becomes exact in the limit 
of a uniform system, the weight factor u must satisfy the normalization condition 

J dx 2 w(xi,x a ;p) = 1 

Requiring that a; must satisfy 

-c l ; (xi,x 2 ;p ) = lim — — — — - (3.5) 
p^po 5p(xi)5p(x 2 ) 



exactly, one finds 

w(xi,x 2 ;p) = 



13 1 c( Xl ,x 2 ;p) + lpAa"(p) 



(3.6) 



2Aa'(p) 

where Aa(p) is the excess free energy per particle and primes on it denote derivatives with 
respect to density. 

The contribution to the free energy due to deformation is given as |I] 

f3AA e {p)=mMp}-^ u {p}) 

= -- / rfxi J dx 2 [pd(xi)p d (x 2 ) - p(xi)p(x 2 )] c(xi,x 2 ;p). (3.7) 

where {3AA u [p] and p(x) are, respectively, the excess free energy and single particle density 
distribution of the ground state of the ordered phase. The subscript d refers to the corre- 
sponding quantities of the deformed state. In writing Eq.(3.7) it is assumed that the direct 
pair correlation functions do not change due to deformation. 

Since the isotropic fluid DPCF is an invariant pairwise function it has an expansion in 
SF frame of the form 
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c(r, Q 1 , n 2 ) = E 51 c (^i> fe, i; "-2; r) C s (Zi, Z 2 , i; m 1 ,m 2 , m) 

hl 2 l rn\mim 

where C g (li, l 2 , l] mi, m 2 , m) are the Clebsch-Gordon coefficients, c(Zi, ^ n ii n 2] r ) the har- 
monic expansion coefficient of the DPCF. f = r/ | r |, is a unit vector along the intermolec- 
ular axis. 

The elastic constants are defined by the second-order term of the expansion of the free 
energy of the deformed state around the free energy of the equilibrium (ground) state in 
the ascending powers of a parameter, which measures the deformation. The first term of 
this expansion is balanced by the equilibrium stresses of the ground state. One defines the 
elastic free energy per unit volume as 

y = ^[AA e [p}+P(V d -V)] 

where Vd is the volume of the deformed sample and P the isotropic pressure. 

For an ordered (liquid crystalline) phase p(x) is expressed in terms of the Fourier series 
and the Wigner rotation matrices |2j as 

p(x) = p(r, fV) = po E E Qimn{G) exp[*G.r] D l mn (n') (3.9a) 

G Imn 

where the expansion coefficients 

21 + 1 r r 

Qimn(G) = j dr J dn'p(r, SI') exp(-zG.r)^ n (n') (3.9b) 

Here {G} are the reciprocal lattice vectors (RLV's) of the positionally ordered structure 
that might be present in the system and po is the mean number density of the system. All 
the angles which appear in Eq.(3.9) refer to a coordinate frame the z-axis of which is along 
the director n. Using the rotational property of the D matrices we rewrite Eq.(3.9a) so 
that all the angles appearing in it refer to space fixed coordinate frame shown in Fig.l. Thus 

p(r,n) = Po J2Y,Y,Qirnn(G)eMiG.r]D l ;jn)Dl n (n) (3.10) 

G Imn q 

In the limit of long wavelength distortion the magnitude of the order parameters are 
assumed to remain unchanged. The changes occur in the direction of the directors making 
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them position dependent and in the RLV's G. The RLV's Gd of the strained structure are 
related to G of the unstrained structure as 

G d = (I + e^G (3.11) 

here e is a strain matrix which governs the change in position. Thus for a deformed state 

p d (r, ft) = PoEEE Qimn(G) exp(tG d .r)D l ; m (h d )D l qn (n) (3.12) 

G Iran q 

where n d (r) indicates the direction of director at r of the deformed state. In the deformed 
state the directors orientations become position dependent. One uses the rotational proper- 
ties of D matrices to re write Eq.(3.11) in terms of the director of a chosen point. Substituting 
Eqs.(3.8)-(3.12) into Eq.(3.7) and after some simplifications we get 

V 2 ° hl 2 l m ' im > 2 m> m\n in2 G (2/ X + 1 ) (2/ 2 + 1 ) 

C^hWMm'X) D l Z x (n) J dr [exp(iG d .r)Z^, m , (A X ) 
- 5 m , m > 2 exp(iG.r)] c(hl 2 l; n x n 2 ] r)Y* m (r) (3.13) 

where n indicates the orientation of the director at the reference point R and Ax the 
angle between the directors at R and R + r. Eq.(3.13) presents a general expression for 
the distortion free energy density for the continuous symmetry broken phases with tilted 
orientational order in the limit of long wavelength distortion. Since for uniaxial phases the 
director is along the z-axis of the SF coordinate frame, D l * lX (n = 0) = 5 m \. When this is 
substituted in Eq.(3.13) it reduces to Eq.(3.9) of I. 

Eq.(3.13) can be used to derive expressions for the contributions arising due to curvature 
in the director orientation and dilation in layer thickness. The term which represents the 
coupling of these two terms of distortions is neglected. Thus 

^ = /. + /. (3-14) 

where 
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(G) Qhm'n 2 (—G) , , 

/i = -«poL L LLL — (9 , t 1V9/ , n — Cgihhi-m^x) 

Z hhl m^m^m' m\ nw 2 G ^1 t lj^(2 t lj 

£>^(n) y drexp(iG.r)[Z^ m ,(A X ) - ^ m ^]^(f) c(hl 2 l; mn 2 ; r) (3.15) 

(G) Ql 2 m' n 2 

^(n) / dr[exp(iG d .r) - exp^G.r)]*^^ YJ^(f) c(ZiZ 2 ^; nin 2 ; r). (3.16) 

While /i includes all contributions arising due to curvature in directors orientations, f 2 
the contribution due to dilation (compression) in layer thickness. 



B. Expressions for the elastic constants associated with directors orientations 

All the elastic constants associated with curvature in director orientations can be derived 
from the expression of f\ given in Eq.(3.15). Using the procedure outlined in II we rewrite 
Eq.(3.15) in terms of the director components. Thus 



1 2v- y y y Qhm' ini (G) Q hm ' 2 n 2 { G) , r 4 r 

1 



dr 



hhlm^m^m 1 m\ nw 2 G (^1 + 1)(2^2 + 1) 

^ | s{ 3) I 2 {(a n + a 12 )(g 2 - g 2 )(r.x)(r.z) + (a 23 + a 21 + a 22 )(g 2 (r.x) 2 + • • ••)} 



+ i I B 2 ] I 2 W«i2 - aii)g 2 (r.y)(r.z) + (a 23 - a 2i - a 22 )(g 2 (r.x) 2 + g 2 (r.y) 2 H )} 



2 

+ 5 I | 2 {^ 2 (g 2 - g 2 )(r.x)(r.y) - ^ 2 (g 2 (r.x) 2 + g 2 (r.y) 2 + g 2 (r.z) 2 ) 

-(an 4 f/ 12 )(//;(r.x)- ----- < ; „, 

"' J 2// 2 „2\ 



- & 2 g^(r.z)(r.x) + ••••} + B ( 2 1) B ( i ) {--(a ll + a 12 )(g 2 (r.x) 2 + g 2 (r.y^ 2 
+ '-(an - a 12 )((g 2 - g 2 )(r.x)(r.y) + ••••) + ^((q'j - ^)(r.y)(r.z) ) 



- ^' 2 (an - a 12 )(g 2 (r.x) 2 + g 2 (r.y) 2 + • • ••)} + 



exp(iG.r)F^(r) c(lil 2 l\ n x n 2 ; r) (3.17) 



where 



a n = -i{(/ 2 + m 2 + l)(/ 2 -m 2 )}2C , 9 (/i/ 2 /;m' 1 ,m / 2 + l,A) 
ai2 = ^{(^2 — + 1)(^2 + m2)}iC ff (ZiZ 2 Z; m' l5 m' 2 — 1, A) 
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a 2 i = + m' 2 + 2)(Z 2 - m' 2 - 1)(Z 2 + m' 2 + 1)(Z 2 - m,)}^^; m' 1? m' 2 + 2, A) 

o 

a 22 = - m 2 + 2)(Z 2 + m 2 - 1)(Z 2 - m 2 + 1)(Z 2 + m 2 )}5C fl (ZiZ 2 Z; m[,m' 2 - 2, A) 

o 

a 23 = -^{(^2 + m 2)(^2 — ra 2 + 1) + (Z 2 - m' 2 )(l 2 + m' 2 + l)}C g (ZiZ 2 Z; m^, A) 

o 

6 2 1 = m' 2 Cg(lil 2 l] m' 1: m 2: A) 

6 22 = m 2 C g (lil 2 l; m' 1: m 2 i A) 

Since the plane wave exp(iG.r) travelling in the ^-direction (because of layered structure 
of smectic phases) is symmetrical about the z-axis and can be expanded as a series of 
Legendre polynomial referred to this axis, 



iGzi2 



^) l \Zl' + l)3v{Gr)P v {cosd) 



(3.18) 



where jv{Gr) are the spherical Bessel functions and 9 is angle between the z-axis and 
intermolecular axis r. 

After performing the angular integration over f we compare the resulting expression with 
Eq.(2.24) and obtain the following expressions for the elastic constants. 



i(2l' + 1) 

m = -pI^k e E E E E(0'' , 2/ V \w 2 , '. Qi^n^Q^i-G) 

hhW m' im ' 2 m\ nm 2 G + L )K Zl 2 + L ) 

'b 2l / 21^+1^ 1/2 



D l * (ft) T nin2 



2l' + l\ 1/2 



C g {2l'l; 000){C g (2l'l; 202)5 m2 - Cg(2l'l;202)5 m2 } 



&22 x 



R (2V + 1\ 1/2 
C ff (2Z'Z; 000) 2 5 m0 + \\A -tzt^t C g (2l'l; 000) 



2 \ 2Z + 1 



(C S (2Z'Z; 202)5 m2 + C g (2l'l; 202)^)} 



(3.19) 



(2/' + 1) 



pb 2 = pIV^Y, E EEE(^ r2/ v , ivoz li^ ^w^^n^-g) 

'b 21 /2Z' + 1^ 1/2 



n<* (n) T nin2 



VQ\2l+l 



C g (2l'l; m){C g (2l'l; 202)5 m2 - C g (2l'l; 202) c^} 



&22 /A /2Z' + l' 

+ -Q-iOu'OmO ~ I "gj 



1\ 1/2 /3/2/' + l\ 1/2 

r J c 9 (2Z'z ; ooo)5 m o - y 2 ( azTlJ 6,9(2/75 000) 



(C S (2Z'Z; 202)5 m2 + C fl (2Z'Z; 202)5 m2 )} 



(3.20) 
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(G)Q hm ' n A-G) 



lll2 U'm [m >mXn in2 G 6 (21, + 1) (2Z 2 + 1) V 

2Z'+ l\ 1/2 



2/ 



C g {2l'l;000) 2 5 m0 



(3.21) 



/3S 13 = -pgV^F £ E E E E(O r , 2/ ~7y7^T~T iT ^-^r^Qi^nxCg) 

hfaW mA G ( 2/ l + X ) ( 2/ 2 + V 2 V6 ^(2/ + 1) 

Qi 2m ' 2 n 2 (-G) D^(n)JJ^,b22C g (2l% 000)[C 9 (2Z'Z; Wl)6 ml + C g (2l'l; 101)* m i] (3.22) 

/L4 = -plV^ E EEE Ui) l ' (2] +fv2/ 1 ^ ^ (Q23 - ° 21 " 

'2Z' + 1\ 1/2 _ _ x9t . /3 /2Z' + 1\ 1/2 



^'<5 m0 - tt^— C S (2Z'Z; 000) 2 5 m0 



2/ + 1; yv ' ' mu v 2 v2/ + i; 

C 9 (2/'/; 000) {C g (2l'l; 202)5 m2 + C g (2l'l; 202)<5 m2 }] (3.23) 



(21' + 1) 

PA 12 = -ftyfc E E E E E(0 ,, f2 , ViV2/ +1) Qi,rn[ m (G)Q l2m/2n2 (-G) 



hhW m\m' 2 mX nm 2 G 

1 /2/' + l\ 1/2 

Dt^KZll^ {-^T) ^ + a i2)G 9 m 000){C g (2l'l; 101)* ml + C g (2l'l; m)S ml } 

1 /2/'4-l\ 1/2 /T /2/' + 1\ 1/2 

+ 3(023 + G 2 i + a 22 ){5 H ,5 m0 - ( 2T+t) ^ (2/,/; 000)2,5m0 + V 2 ( 2ZTT ) °° 0) 

(C ff (2Z'Z; 202)5 m2 + C 9 (2Z'Z; 202)5 m2 )}] (3.24) 

/f2Z' + 1) 

pa 21 = -piv^ e E E E E(0 / ' 7 ^ w ^^ im>1 (G)g i2m; „ 2 (-G) 



l\l 2 ll' m^m'z mX mn 2 G 



1 /2/' + l xl/2 



JhW'[~^ {-2TT1) ( ° 12 ~ ail)C ^ 2l% 000 ){G g (2l'l- Wl)5 ml - C 9 (2l'l- m)5 ml } 

+ 3(023 - a 2 i - a^HMmO - ( 2ZTT ) C9(2/7; 00 °) 2<5m0 " V 2 ( 2ZTT ) °° 0) 
(C 9 (2Z'Z; 202)5 m2 + C g (2Z'Z; 202)5 m2 )}] (3.25) 

PC, = ^ E E E E ECO^ ToT^^^ Qi^n^^Qh^C-G) 

2 hhll' m' x m! 2 mX n\n 2 G (2ll + L)(2l 2 + L) 



1 /2/' + l\ 1/2 

D l * x (n)J%%,[^ ( 2ITtJ (an " a ^ C ^ 2l ' l i 000){C 9 (2I'I; 202)5 m2 - C g (2l% 202)5 m2 } 

1 1 , (2V + l\ ^ 2 

- g (an + 012) + m' 2 (a u - a 12 )}{5 w 5 m Q - ( - - - J C g {2l'l; 000) 2 5 m0 

R /or + 1\ 1/2 

+ J- UtTT C , 9 (2/7;000)(C , ,(2/7;202)5 m2 + ^(2/7;202)5 m 2)}] (3.26) 



f3C 2 = - P0 ^C~ J2 EE , 9/ _/ n W9/ i n^im' 1 n 1 ( G ) ( 5« 2 m;n 2 (- G ') 



1 /2/'+l\ 1/2 

Dtx(n)Jim^ {-n-lTT) ^ ~ "nWgWh 000){C 3 (2/7; 202)5 m2 - C S (2Z7; 202)5 m2 } 
+ g ( a n + a i2) + "i 2 ( a n - a 12 )}{5w5 m Q - ( 2l + 1 J C g {2l'l; 000) 2 5 m0 

~\Z|(iTl) 7 ^( 2 ^ 000 )(^( 2 ^^^^ 

1 /2/' + l\ 1/2 

- -^m 2 C fl (2i7;m' 1 m 2 A) ( J C g (2l%000){C g (2l'l; 101)5 ml - C 9 (2/7; 101)5 m i}] (3.27) 



where 

J?ML = Jdrr 4 j L (Gr) c(hl 2 l; mn 2 ; r) 

are the structural parameters. A line under a numerical subscript denotes a negative quan- 
tity. Eqs.(3.18)-(3.26) give general expressions for the curvature elastic constants of the Sm 
C and other phases of this class in the long wavelength limit. From these results one can 
obtain expressions for the elastic constants of a given phase by using the appropriate order 
parameters corresponding to the phase. 

The structures of smectic phases are characterized by three different class of order pa- 
rameters; (i) orient at ional, (ii)positional, and (iii) mixed ||. In general there can be upto 
{21 + l) 2 orientational and so many mixed order parameters of rank I for a given value of 
the reciprocal lattice vector G. One uses the symmetries of the phase and of the constituent 
molecule to reduce this number. Since SmC phase is biaxial with point symmetry C 2 h, we 
choose the following orientational order parameters at / = 2 level (molecules are assumed to 
be cylindrically symmetric). 
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A = ^Q 2 O = (^(COS0)) 


^1 = ^Q2i = -7Q21 = (7^) (sin 9 cos (9 cos 0) (3.28a) 



5 - 5" V2 



w 2 = -Q22 = -Q22 = (-) 2 (sin 2 6 cos 20) (3.28b) 

(3.28c) 

While P 2 measure the degree of alignment of molecular axis along the director n which makes 
a non zero angle with the smectic plane normal, iO\ and u 2 measure the degree of alignment 
along the director c and the biaxial ordering. Angular bracket indicates the average over 
the orientational distribution . 

For positional ordering one may choose as in the case of Sm A phase, one order parameter 
corresponding G = G z = ^ where d is the layer spacing in equilibrium. Thus 

H = Qvo(G) = (cos2tt^) (3.29) 

a 

The coupling between the positional and orientational ordering may be described by the 
order parameters Ti, r 2 and r 3 , where 

1 2ttz 

ti = 7 Q2o(G) = (cos— P 2 (cos#)> (3.30a) 
5 1 a 

r 2 = -Q 2 i(G9 = ~-Q2i(G) = (-)i(cos^sinflcosflcos0) (3.30b) 
5 5 la 

r 3 = \Qn{G) = \q 22 {G) = (^(cos^sin 2 #cos20> (3.30c) 
o a 

Using these order parameters we give below explicit expressions for the curvature elastic 
constants for Sm C phase consisting of cylindrically symmetric molecules; 

B 1 = M 1 sin 2 ip + M 2 cos0> + M 3 (3.31a) 

B 2 = M 4 sin 2 i) - M 2 cos i/j + M 5 (3.31b) 

B 3 = M 6 sin 2 if> + M 7 (3.31c) 

B 13 = (3.31d) 

A = M 8 sin 2 ^ + M 9 (3.31e) 
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A 12 = Mio sin 2 ip + Mn (3.31f ) 

A 21 = Mi 2 sin 2 V + Mi 3 (3.31g) 

Ci = M 14 sin 2 ^ + M 15 (l + cos^) 2 + M 16 (l - cos^) 2 + M 17 sinV> + M 18 sin ip cos ^ + M i9 (3.31h) 

C 2 = M 20 sin 2 ip + M 2 i(l + cos^) 2 + M 22 (l - cos^) 2 - M 17 sin^ + M 23 sin?/' cos ^ + M 24 (3.31i) 

Here ip is the tilt angle. The directors n and c of a uniform phase are assumed to be 
in the xz plane. The expressions for the constants Mi — M 2 4 which involve the order and 
structural parameters and are, therefore, temperature and density dependent are given in 
Appendix A. In writing these expressions we have taken into account the fact that for a 
system of rigid axially symmetric molecules n\ = n 2 = and 

1 



c(hl 2 l] 00; r) = — yj(2h + l)(2/ 2 + l)c hhl (r) (3.32) 
The structural parameters are connected with harmonic coefficients ci x i 2 i as 



J hl2l (0) = -?L [°° drr 4 c hhl (r) (3.33) 



2 

4, 



Po ro ° 



and 

J hlaW = / drr A ]v {Gr)c h Ur) (3.34) 

V47T J 

Our results given in Eqs.(3.31) show that for cylindrically symmetric molecules = 0. 
This is, however, not true in general. For molecules with broken cylindrical symmetry Bi 3 
will be non-zero. Its magnitude, however, appears to be small compared to Bi,B 2 and B%. 
Madhusudana and coworkers |§ have found ~ 0.3 



C. Expression for the Elastic constant B 

As mentioned above the expression for the elastic constant B is derived from the expres- 
sion of f 2 given by Eq.(3.16). Since Sm C as well as Sm C* has one dimensional positional 
ordering in addition to orientational ordering we have 
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r 1 Qlirn\ni(G)Qi 2m ' n2 ( y G) , . 

/2 = -^po L L L L L — r9/ t 1V9/ Vt; — ™i™ a) 

2 G iit 2 l m^m' mA nma 1^1 + 1 J 1^2 + 1 J 

D l * nX (h) J dr [exp(iG d 2i 2 ) -exp(iGzi2)]<J ro 'mic(ZiZ2Z;nin2;r)YJ^(f) (3.35) 

where G = ^ and G rf = ^p; d e and d are the interlayer spacing of the distorted 
and undistorted Sm C (and Sm C*) phase. Using eqn.(3.11) we expand above equation in 
ascending power of dilation parameter e = ^ — 1. The term associated with |e 2 defines the 
bulk elastic constant B. Comparing Eq.(3.35) and with Eq.(2.4) we get, 

B = yPo LL L LL — (9 , T 1V9/ , n — Cgihhl-m.mX) 

1 G hhlm'^m' mX nma ^ + W'2 + lj 

£>L(A)G 2 / rfr 4 exp(iGz)c(ZiZ 2 Z; mn 2 ; r)Y^(r) (3.36) 
Using Eq.(3.18) we have 

_ kTpQ ^J~^ y^ y-^ y-^ Qhm'^rn (G)Q; 2m ^ W2 (~G) 



6 C tWI'SraX A (2Zi + l)(2Z 2 + l) 



G 2 (2Z' + l) 1 / 2 D^(n) C^Z.Z; m! x m!\)J™ 



/2Z' + 1\ 1/2 

5 »' + 2 (^rj Cfl(2n;000) s 



(3.37) 



As shown in I, the term proportional to e in the free energy expansion gives the condition 
for the interlayer spacing. In terms of the order parameters given by Eqs.(3.29)-3.30), 
Eqs.(3.37) can be written as 

(3B = ^ sin 2 ip + L 2 sin -ip cos ip + L 3 (3.38) 

The expressions for Li, L 2 and L 3 for a system consisting of cylindrically symmetric 
molecules are given in Appendix B. 



IV. DISCUSSIONS 

In order to estimate the values of elastic constants given above we need the values of the 
order parameters, generalised spherical harmonic coefficients of the direct pair correlation 
functions of an effective fluid as a function of temperature and density and the information 
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about the constituent molecules, viz. electric multipole moments, geometry of the repulsive 
core, length-to-width ratio, etc. as input parameter. 

Since in the limit of long wavelength distortions the magnitude of the order parameters 
are assumed to remain unchanged, the value of the order parameters at a given temperature 
and the density may be taken either determined experimentally or calculated from a theory. 
The c-harmonics for a given system can in principle, be found either by solving the Ornstein- 
Zernike equation with suitable closure relations || or by adopting a perturbative scheme 
which is based on the fact that the fluid structure at high densities is primarily controlled by 
the repulsive part of the interactions. However, such calculations for non-axial molecules are 
very complicated and may need enormous computational efforts to generate reliable data 
for c harmonics ||. 

The difficulty which arises in applying the theory to real systems is related to the potential 
energy of interaction between mesogenic molecules. The mesogenic molecules are large and 
have groups of atoms with their own local features. One way to construct the potential 
energy of interaction between two such molecules is to sum the interatomic or site-site 
potential between atoms or between interaction sites. However, for mesogenic molecules 
there are too many terms in this sum to be practical. Moreover, the dependence of interaction 
on molecular orientations in this expression is implicit so that it is difficult to use it in the 
calculation of angular orientation which give rise to liquid crystals. 

In another and more convenient approach one uses rigid molecular approximation in 
which it is assumed that the intermolecular potential energy depends upon only on the po- 
sition of the centre of mass and on their orientations. This kind of approach neglects the 
flexibility of molecular structure which plays important role in stability of many liquid crys- 
talline phases. In view of various complexities in the intermolecular interaction one is often 
forced to use a phenomenological description, either as a straight forward model unrelated to 
any particular physical system, or as a basis for describing by adjustable parameters fitted 
to experimental data for interaction between two molecules. Most commonly used models 
are hard-ellipsoids of revolution, hard sphero- cylinders, cut-sphere ]10| and Gay-Berne [11 
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However, non of these are known to show the existence of Sm C phase. 

The computer simulation study [HJ shows that the GB model exhibits Sm A and 
Sm B phases. The GB potential contains four independent parameters which control the 
anisotropy in the attractive and repulsive interactions and can be written as 

O"0 ^ 2 r °0 



u(l, 2) = 4e(r , n 1? O a ) 



{- 



12 

-} "{- 



"I 



(4.1] 



r - a(r, f2 1( fi 2 ) + 0"o r - o"( r , ^2) + Co' 

Here e(f, f2i,f2 2 ) and a(r , f2i, f2 2 ) are angle dependent strength and range parameters, 
respectively, and are defined as 

,(r.e0 2 + (r.e 2 ) 2 - 2 X '(r.e 1 )(r.e 2 )(e 1 .e 2 ) 1 2 



6(r,O 1 ,O 2 ) = e [l-x 2 (ei.e 2 ) 2 ]- 1/2 - \' 

(r.ei) 2 + (r.e 2 ) 2 - 2x(r.ei)(f.e a )(ex.e a ) 



cr(r, Oi, Q 2 ) = cr 



l-X ,2 (ei-e 2 ) 2 

-1/2 



i-x : 



(4.2) 



1 - X 2 (e"i.e 2 ) 2 

where ei and e 2 are the unit vectors along the symmetry axes of two interacting molecules, 
eo and <jq are parameters which provide a measure of the attractive interactions and the 
width of the molecule. The anisotropy parameter x an d x' are defined as 



X 



xl-l 1 k' 1/2 - 1 

and x ~ 



xl + l ^ A fc'Vs + l 
where xq = 2a /2b is the length (major axis) to breadth (minor axis) ratio and k' is the ratio 
of the potential well depth for the side by side and end to end configurations. We have 
taken here Xq = 3.0 and k' = 5. The harmonics were generated by solving the Ornstein- 
Zernike equation using the Percus-Yevick closure relation |T3| . Using these c-harmonics we 
have calculated the values of the structural parameters at ksT/e^ = 0.8 and the packing 
fraction rj = (|)xoPo°"o = 0-49. The values of few structural parameters for G = and 
G = — = 2.0944 so obtained are as follows 

J 220 (0) = 2.23174 Jqooo = -2.78902 

J 222 (0) = 0.33 3 24 J0002 = -4.61192 

J 224 (0) = -4.02296 J 202 o = 0.07539 

J 242 (0) = 0.30 951 J 202 2 = -2.05 8 25 
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J 2 44(0) = 


-1.06209 


J2024 


= -0.56503 


^440(0) 


= 0.49893 


J220Q 


= 0.13322 


^442(0) = 


-0.09114 


J2202 


= 0.6818 


^444(0) 


= 0.13077 


J2220 


= -0.49332 








= -0.00627 






J222A 


= 0.17776 






J 22m 


= 0.54022 






J 2^20 


= -0.25849 



Since the values of the order parameters are not known, we assume them to be equal 
with value 0.5 to estimate the relative contributions made by different terms in Eqs. (3.31) 
and (3.38). It is found that (i) dominant contributions for all elastic constants come from 
angle independent term, (ii) coefficients of linear terms i.e. terms involving simp or cosij) 
are small compared to terms involving sin 2 -?/>and (iii) the relative contributions to A, A± 2 , 
A21 made by angle dependent terms are small compared to their relative contributions to 
By, B 2 and B 3 . The values of all these constants except Bi 3 are found to be of the order 
of 10~ 6 dynes. As has already been pointed out that for linear rigid molecule B 13 is zero. 
For molecules with broken axial symmetry _B 13 is non-zero but its value is less compared to 
other elastic constants. The elastic constant B is found to be of same order as in Sm A. 
However, as tilt angle increases the value of B decreases. The two angle dependent terms 
in Eqs. (3.38) is found to cancel each other, and therefore, the major contribution is due to 
the angle independent term. 
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APPENDIX A 



In this appendix we give expilicit expressions for the coefficients M 1 to M 2 4 in units of 
ksT/aQ (see Eqs.(3.31)). As is obivious from Eqs. (3.19)-(3.27), they involve order and 
structural parameters of all orders. We list only first few non-vanishing terms. 



Mi = -u\ 



M 2 = to 2 



M 3 = lo\ 



M 



1 5 



J 222 (0) 



2 V 14 

^222(0) 



LUn 



,10 

0WyJ 222 (0) 



14 

:i /i 



U/ — J 222 (0) 



3 V 14 



- P 2 UJ 2 

(-2J 2 2 2 (0) + Vl4J 220 (0)) 



+ P2UJ2 



+ 



4 15 



J 222 (0) 



+ 



+ u 2 2 



4^(J 22 o(0)- 1/^222(0)) 



+ 



+P 2 u; 2 



hi— J 222 (0) 



+ 



4 = -T1T3 



30 10 



M 5 = uj\ 



7 V 7 

1 /5 



J 5 



2222 



+ 



3 V 14 



(4J 222 (0) + V14J22 (0) 



LUn 



4^(J 2 2o(0) - 



^222(0)) 



-P 2 CJ 2 



U/— J 222 (0) 



+ 



M fi = 



U/yJ222(0) 



M 7 = -u; 2 



M s = -P, 2 



M 9 = uj\ 



1 no 
3vY 

To 
7 

2 /5 



J 222 (0) 
^222(0) 



+ ^2^2 

/7 



4 15 



3V 7 

^220(0)) 



J 222 (o) 



+ 



4^(J 2 20(0)+ 2J^J 2 22(0) 



+ 



+ 



a/5 

, v L4 J 2 2 2 (0) + ^-J 22 o(0) 



-P 2 

r 2 



3^ 



(J 220 (0) 



2 /2 
3V7 



J 222 (0)) 



+ 



M10 = 



-P 2 

-*2 



10 T 

y^222(0) 



+ 



Mn = ^1 



M12 = 



^^(-14J 222 (0) + J 220 (0)) 



+ Po 



yU 20 (0) 



+ 



15 10 

YVT 2222 



+ 
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1 



M 13 = uj\ 

M15 = h^i 

M 16 = P 2 CU! 



21 V 14 

35 / 3 



(28J 222 (0) + 7^14/220(0)) 

^222(0) 



"V5 



^220(0) 



+ 



8 V 35 
25 



16-v/l05 
5 



J 2 2 2 (0) 

•Wo) 



+ 
+ 



M 17 = -P 2 2 



M 18 = P| 



1 / 5 



+ 



Wn Ja22(0) 



Mig = P 2 ^i 
^20 = -A^l 
M 2 l = -P 2 Wl 
M 22 = —P2W1 

M 23 = 7? 
M 24 = P2W1 



|JA(7^(0) + f,/l.W0)) 



+ 



5 / 3 . . 

sV^ 222( } 

S\/I j222(o) 

48V35 J - (0) 



+ 



+ 



15 5 
28 V 14 32222 



+ 



4V|( 7 ^2(0) + ^^J220(0)) 



+ 
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APPENDIX B 



We give explicit expressions for the coefficients L x — L 3 involving first few terms of order 



and structural parameters. 


6 ^ 

u G 


- 3)GVl/022 + + ^^^7222 + • 


6 ^ 

u G 


\/QG 2 UToInoo + A -(3V2 - \/3)G 2 TiTiIoio + 

V v ' ^ ' f-v > 2-*- \)22 * W y\ u V ^ \ tj jv^t 1 ^ 1 2 1 - 222 * 


1 

u G 


/in 

Ai^/ooo + 2GW022 + v^GV^o ~ JyGV? 


where 




-^ooo — ^0000 " 


- 2 J0002 


-^022 — 2Jo220 


55 36 

Y '0222 + "^-<A)224 


-^220 — J r 2200 _ 


- 2 J r 2202 


^222 = 2J2220 


55 36 

y-^2222 + "y-^2224 
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FIGURES 

FIG. 1. Molecular arrangements in Sm C, the directors n, c and the space-fixed coordinate 
frame 
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This figure "ysl.gif" is available in "gif" format from: 
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